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A NEW LOOK AT DECOMPOSITION OF TURBULENCE FORCING FIELD AND 
THE STRUCTURAL RESPONSE 


Y. K. Lin and S. Maekawa 

University of Illinois at Urbana-Chainpaign, Urbana, Illinois 

Abstract 

Measured cross-spectrum of a turbulence field usually shows some 
decay in the statistical correlation in addition to convection at a 
characteristic velocity. Under such a random excitation the computation 
of structural response statistics becomes much more tedious than that 
which would be the case if the turbulence were convected without decay; i.e. 

I 

convected as a frozen-pattern. It is shown in this paper that a decaying 
turbulence can be decomposed into frozen-pattern components thus permitting 
a simpler way to calculate the structural response. The procedure so de- 
vised also provides a relationship whereby the measured input spectra can 
be incorporated. For illustration the theory is applied to an infinite 
beam which is backed on one side by a fluid-filled cavity and is exposed 
on the other side by the turbulence excitation. The effect of the free 
stream velocity is also taken into consideration. 


A NBW LOOK AT DBCOMl’OSITION OF TUUBULENCH I’ORCING FinU) AND 
TMG STRUCTURAL RnSFONSE 


Y. K. Lin* and S. Maekawa// 

University of Illinois at Urbana-Chainpaign, Urbana, Illinois 
Introduction 

From the standpoint of structural response calculation the simplest 
mathematical model for an atmospheric or boundary-layer turbulence is one 
that is statistically homogeneous in space and is convected in a given 
direction as a frozen pattern. The second part of the assumption is known 
as Taylor's hypothesis which often results in tremendous computational 
savings. In some cases it may be the key assumption making tlie problems 
solvable. 

However, experimental measurements of real turbulences invariably show 
that spatial decays do exist in the cross-spectra or cross correlation func- 
tions. Such decays are indicative of the change in turbulence patterns as 
they move down-stream. Tlierefore, results obtained from a frozen-pattern 
analysis are just crude estimates which stand to be improved when a better 
method becomes available. 

The primary objective of this paper is to show that a decaying turbulence 
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can bo constructed by superposing infinitely many frozen-pattern components 
with random amplitudes and convected at different velocities. Then the 
structural response can also be similarly superposed with each component 
corresponding to one frozen-pattern component in the forcing field. 

To focus attention on the essentials, the discussion in this paper 
will be restricted to one-dimensional space coordinate. We shall begin by 
reviewing some basic relations for the frozen-pattern case as building blocks 
for the later superposition. Then the proposed superposition will be devel- 
oped. Finally, as an example, the theory will be applied to an infinite 
beam under the excitation of a supersonic boundary-layer turbulence. ' . 


Frozen-Pattern Turbulence 

An example of one-dimension'al structure exposed to the excitation of a 
turbulent pressux'e field is depicted in Fig. 1. The x coordinate frame is 


stationary with respect to the undeformed structure, and it will be referred 

to as the fixed frame in the sequel. If the pressure is truly of a frozen 

tyiJe and is convected at a constant velocity in tlie positive x-direction, 

then it is a random function of x - U t. Such a random function can be cx- 

c 

pressed as a Fourier-Stieltj os integral as follows: 

pcx-u^t) = i <ink) 


Cl) 


whei'e the frequency w and wave number k are related to the convection speed 
U as w/k = U . It is known from the random process tlxeory that 

O 0 


E (dFCkp dF*(k2) } = SpCkp 6(k^ - k2) dk^ dk2 


C2) 


where n{ } represents tlio ensemble average, an asterisk denotes the complex 


conjugate, and S^^Ck) is the wave-number spectrum in a coordinate frame moving 
at the velocity (referred to as the moving frame in the sequel). 

The cross-correlation function E {p(Xj^ - U^t^) p (^^2 “ ^c^2^^ pres- 

sure, referi'ed to the fixed frame, can be calculated simply by use of Eqs. (1) 
and (2). This function, denoted by Up, depends only on ^ - U^T where ^ “ Xj^ - 

and T = tj^ - and it is related to the moving-frame wave-number spectrum 


S (k) as follows: 
P ^ 


RpC? - UJ) = 


eik^V “ SpCk) dk 


(3) 


If a Riemann-Fourier transform is taken of Eq. (3) we obtain the fixed- frame 
frequency cross-spectrum of p: 




R <5 - u T) dT 

p C 


P V 


C4) 


Equation (4) shows that the fixed-frame frequency cross-spectrum of a frozen- 
pattern turbulence has a special form where ^ appears only in the imaginary 
exponent. This equation also provides a simple formula to convert from Sp(kJ 
to «J)pCC, w) . 

Conversely, to convert from w) to S (k) ; 

P P 


Sp(k) = (‘'cl ‘V 


(5) 


Evaluated St ^ = 0 the cross-spectrum 4*pC5j w) reduces, of course, to the 

m 

usual spectrum. 

We emphasize that Eqs. (4) and (5) are valid only if the turbulence is 


strictly of a frozen pattern, and is convccted at speed U^. 

Equation (1) suggests that structural response to a frozen-pattern 
turbulence can be constructed from a fundamental solution where the excita- 
tion is just a convected sinusoidial pattern of unit amplitude. Tluis, let 
HCxj k) oxpCitot) be the steady-state solution for 

{ HCx, k) Qxp(iWt) } « exp [iCwt - kx)] C6) 

where, symbolically, 0^ represents a linear operator in x and t, pertaining 
to the dynamic problem at hand. Of course, tliis solution must satisfy all 
the necessary boundary conditions. Tlien tlie solution to 

oC^ wCx, t) } = p(x - U^t) • (7) 


after reaching stochastic stationarity, may be expressed as 


w(x, t) = 


H(x, k) expCiWt) dFCk) 


= i: 


H(x, kj e.xpCiU^kt) dF(k) 


( 8 ) 


It follows that the cross-correlation function of the structural response 


is 


E{ w(x , tj wCx,, t )} 

00 

H(x^, kp H*(x 2, k2) e ^^c^^l^l " S^Ck^ 6Ckj - k^) dk^ dk2 


HCXj^, k) li*Cx 2 . k) e ^*^^1 " S fk) dk 


(9) 


As expected, this correlation function is dependent only on tj^ - t 2 « If it is 
desired to calculate this correlation function in the frequency domain, we may 
substitute Eq. CS) into Eq. (9) and change U^k to w. 
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U{ w(Xj^, tj^) wCx^, t^)} 

ll(Xj., w/U^) H*CX2. w/iy “ ^2^<I<p(0, w) clw (10) 

J .00 

In terms of the input and output spectra the relations are eJ^tremely simple 
and illuminative; they are: 

in the wave number domain: 

X2; k) = IICXj^, k) H*(x 2, k) Sp(k) (11) 

in the frequency domain: 

x^i to) = H(x^, (O/U^) H*(x 2 , to/U^) $p(0, to) (12) 

When Xj^ = x^ these formulas reduce to those for the usual spectra, and they 
have the same form as the well-known result for a single degree of freedom 
system in the random vibration theory. The simplicity is a direct consequence 
of the frozen-pattern assumption. 


Decaying Turbulence 

Measured frequency cross-spectra for real turbulences with respect to a 
fixed frame of reference have the general form of 

$ (C, to) = ¥„(0, to) (O exp(-ilo?/U ) (13) 

where is an even non-negative definite function of having an absolute 
maximum equal to one at ^ = 0 and approaching to zero at large absolute values 
of This form is sometimes attributed to Corcos^. A number of researchers 
have I'eported curve- fitted results for (0> t*J) snd (^) . ' For representative 

Jr 

2 3 

works we cite the papers by Bull , Willmartli and Wooldridge , and Maestrello* et.al. 
Implicit in Eq. (13) is tliat a real turbulence is not a frozen one. 


To obtain a theoretical spectrum consistent Kitli Uq. (i3) we propose 
tlie following representation for a general turbulence pressure 


p(x, t) = 


W 


(143 


P(x - ut) clG(u3 

.CO 

Equation (14) implies that p(x, t) is a superposition of infinitely many 
frozen-pattern components, each having a random amplitude dG(u) and a con- 
vection velocity u. Such velocities con assume either positive or negative 
values. Of course, oacii frozen-pattern component can, again, be decomposed 
into sinusoids. Thus, 


P(x, t) = 


^^iCupt - px) 


(15) 


and its fixed- frame correlation is 
E{ p(x^, tp p(x2> t^)} 

E{ dFCP^, up dF*(32^ up dG(up dG*(up} (16) 

In order tliat this correlation function may depend only on ? = x^ - x^ and 
T = t^ - t-, which we sliall assume to be true, the ensemble average under the 
integral sign in Eq. (16) must have the form 

E{ dFCP^, up dF*C32> ^2^ dG*(up} 

= SpCpj, up 6(3^ - 3p iS(u^ - up d3^ d32 du^ du2 

Substitution of (17) into (16) results in 


(17) 


R,(?, T) = 


p(3uT - 35) u) 


(18) 


We nov; apply a Fourier transformation to obtain the fixed-frame frequency 


spectrum 
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1 -itOC/u c; 

I —1 0 £> 


OJ I 

p u du 

^ U 
^ / 


Cleiirly Eq. (19) is a generalization of Eq. (4). 

To compare Eqs. (19) and 13), the latter is Fourier-transformed 


to yield 


where 


4'p(w. 0 <15 

= ¥p(0, w) T (a - w/U^) 


#co 

T Cv) = -|^ (?) e^^'' d? 

J wCO 


Therefore, 


$ (?, w) = $ (0, w) T (a - w/U ) e ^ da 


Letting a = oi/u we obtain 


<5* (?, w) = 


CO, uof I4» c^-f) 

■ -00 u c 


Then equating and we find a formula to compute S (w/u, u) 


as follows; 


SpCio/u, u) = [w/u| T (w/u - co/Uj *p(0> w) 


P J_oo 


Hie frequency cross-spectrum for the structural response can be ob- 
tained by a similar superposition. Thus by a generalization of Oq. (12), 


f 




f I 


• S- 


^2' H*(^2’ Sp(w/u, u) du 

Ox', lotting k a w/u. 




\i^y ^2> = 




Now since 


SpCk. w/k) =lk|?p(0, w) H' ( k - w/ty 
Xv'o obtain a very simple result 

X 2 ; to) = ? (0, to) IlCXj, k) I 1 *(X 2 , k) V (k 


to/U^) dk 


(25) 


(26) 


(27) 


An Example 

As an example, the theory will now bo applied to an infinite beam shown 
in Fig. 1. The beam is backed on the lower side by a space of depth d which 

I 

is filled with an initially quiescent fluid of density p 2 tmd sound speed a 2 > 

On the upper side tlie beam is exposed to the excitation of a supersonic 
boundary- layer turbulent pressure p. The fluid on the upper side of the beam 
which carries the turbulence lias a free stream velocity U^, density and 
sound speed a^. 

As the beam responds to the excitation its motion will generate additional 
pressures in the fluid media on the upper and lower sides. Denoting these 
generated pressures by p^^ and P 2 > respectively, the governing equation of the 
beam motion is given by 

El ,n = p + (p. - p,) (28) 

ay at^ ^ 

For the purpose of determining tlie "wave-number response function" 

H(x, k) the turbulent pressure p should be replaced by exp[ i(wt - kx)] and 


i 


) 


>9 


tlic structui’al response w equated by ll(x, k) cxpCiwt) c A(k) cxpC-ikx) 
expCiwt). Furthermore, we shall make tlio usual approximation tliat pj^ can 
be calculated without regard to the presence of the turbulence. Then p^ 
is governed by the equation 


,9 .. 9.2 2 , 9 ^ 9 ^ . 

*- 9 t ^«> 9 x ^ ^1 " ^1 ^ g ^2 ■*“ g ^ 2 ^ ^^1 " ° 


(29) 


and subject to the conditions that p^^ can propagate only in the positive z 
direction, and that 

s 

The solution for p^, when evaluated at z=0, is known to be 


3Pi 


^9z •' 


z=0 


-ikx iwt 


kCw/k - UJ" 

CPi),-0 2 ® ° 

^ ^ i Wk^uX - F 


(30) 


The pressure generated on the lower side of the beam is governed by 


the equation 

9^p. - .2 r,2 

^2 „ 2 ,0 

2 " ^9 C 2 ^ 9) P9 " 0 

9t^ 9x^^ 9z^ ■ 

and subject to the conditions 
3P2 

— — =0 at z = -d 


(31) 


and 


3p. 

9z 


w at z=0 


The solution for P 2 > when evaluated at z=0, is given by 

^ ^ .2 cotyd , ^ “ikx iwt 

^Vz-0 " ^' 2 ^' -y-A(k) e e- 


( 32 ) 
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whox'c 


liquations (30) and (32) can now bo substituted into liq. (28) to find 

A(k), recalling that p must bo replaced by exp[ i(wt - kx)] and w by A(k) 
••ikx 

c ' e . The result may be expressed as 
A(k) = 1/Q 


where 


Q = BIk'' - mw^-h ipjOj - UJ ^ 

[cw/k - uy - 

+ (cotyd)/Y 


(33) 


71ius , 

H(x, k) = exp(-ikx)/Q (34) 

We arc now ready to compute tlie structural response spectrum due to the' 
excitation of a convected decaying turbulence. For the input spectrum of 

4 

the turbulent pressure, Eq. (13), we shall uso a form proposed by Mac-strollo , 
for which 

t|^(C) = oxp( “*^ 3 

where 6 is the boundary- layer thickness and a is an experimentally detei'mined 
quantity. Corresponding to this i|/ function we have 


'J' (k - oj/U ) = {Tra6 [Ca6)“^ + (k - 

c c 

Having determined H(x, k) and T (k - w/U^) the cross-spectrum of the structural 
response may bo computed using Eq. (27). T)ie integration -over k must be 
carried out, however, on a digital computer. 

Figure 2 shows the computed results for the frequency spectrum (i.c., 
when Xj = x^) of the structural response using the following physical data: 


-II- 


h 



proportios of the boom 

4 2 

El (bending rigidity) n 3.945 x 10 N-m 
in (mass per unit length) e 9.746 kg/m 
properties of the surrounding fluid media 

3 

» P 2 » p (density) = 0.11015 kg/m 
^1 " ^2 ~ ^ (speed of sound) 261. G m/scc 

(free stream velocity on upper side of beam) « 575.6 m/soc 
U (convection velocity of the turbulence) = 0.75 
d(cavity depth) = 0.1178 ra 
properties of t!ie turbulent pressure ^ 

I s! 

rjj(0 = decay factor = exp ( - ) 


^’p(OjW) = spectral density 


2U, 


r A 


M/UJ 


00 n=l 


n 


6 (boundary- layer thickness) = 0.279 m 


and experimentally determined constants 


0 = 3 

= 4.4 X 10"^ 
A2 = 7.5 X 10"^ 
A^ = - 9.3 X 10"^ 

A, = - 2.5 X 10“^ 
4 


« 5.78 X 10 ^ 
K 2 = 2.43 X 10“^ 

Kj = 1.12 


= 11.57 


Concludin}’ Roiiiai’ks 


The theory tlcvoloppd liorein is applicable to any turbulence forcing 
field which has a cross-spectral density of the form of liq. (13). It is 
particularly useful in dealing with boundary- layer turbulence for which 
the decay in correlation is much more significant than that of the 
atmospheric turbulence as far as structural response is concerned. For 
this reason we have chosen an example to illustrate tlie application of 
the theory whicli includes the effect of a cavity and tlie effect of the 
free stream velocity. These are main features in the problem of fuselage 
panel vibration under the excitation of boundary-layer turbulence. The 
infinite unsupported beam is perhaps the simplest structural model possible 
which still allows consideration of these features. The advantage of a 
simple model is to avoid the burdensome iiiathematical details and concentrate 
on basic principles. . 

A better representation of fuselage panels can be obtained by adding 
eve'nly spaced elastic supports to the infinite beajn. The model then be- 
comes a periodic structure; i.e. , a structure which is composed of identical 
sub-units and for which analytical studies have been carried out exten- 
sively^’^'^^^’^*^'^^’ The elastic, supports give rise to multiple reflections 

and the solution beco-ii.'is considerable more complicated. Details of this 

12 

solution will be reported in another paper 

Further extensions to the two-dimensional case are obvious. With 
small modifications the solution for the periodic beam can be changed to suit 
the case of a row of panels simply supported along two parallel edges. In 
the sense of Levy’s series solution to plate problems involving two simply 


' - 13 - 

sujDportcd parallcil edges, tlio solution is still mathematically exact. 
However, if tlic structural model has more than one row of panels then 
only approximate solutions are possible since the spatial variables now 
cannot be separated. Although now concepts are not required in treating 
the two dimensional problems the machine computation time can increase 
astronautically and such studios may best be cai’ried out in the industry. 
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figure 1. All infinite beam under the excitation of boundaiy- layer turbulence, 
figure 2, Specti'al density of structural response, 
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p (turbulent pressure) 





